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Abstract 


A procedure to recover explicitly self-adjoint matrix Dirac sys¬ 
tems on semi-axis (with both discrete and continuous components of 
spectrum) from rational Weyl functions is considered. Its stability is 
proved. GBDT version of Backlund-Darboux transformation and var¬ 
ious important results on Riccati equations are used for this purpose. 
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1 Introduction 


Self-adjoint Dirac system has the form 



( 1 . 1 ) 


where 


/mi 0 

0 -Im, 


V 


0 V 
V* 0 


rui -|- m 2 ='■ m, (1-2) 
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Jmj. is the rrik x identity matrix and v{x) is an mi x m 2 matrix function. 
We assume that the potential v is locally summable (i.e., summable on all 
the hnite intervals [0,/]). 

Inverse spectral problem for general-type self-adjoint Dirac system, and 
closely related problem to recover Dirac system from its Weyl-Titchmarsh 
(Weyl) function, had been actively studied since 1950s [28, 32] and many 
interesting results were published last years (see, e.g., [2,3,6,7,10,23,35,43] 
and various references therein). Further, speaking about inverse spectral 
problems, we mean (in particular) the recovery of systems from their Weyl 
functions. Inverse spectral problems (in the mentioned above sence) are also 
solved [12,14,37,39,43] for general-type skew-self-adjoint and discrete Dirac 
systems. 

Procedures to solve these inverse problems are nonlinear and usually un¬ 
stable. However, stability plays an essential role in theory and applications, 
and several special cases where such stability could be proved are important. 
In particular, we could mention the paper [34] on the evolution Schrodinger 
equation and the paper [23], where stability was proved for a class of scalar 
(mi = m 2 = 1) Dirac systems on interval with discrete and d-separated spec¬ 
tral data. Here, we consider the case of explicit solutions of inverse problems 
(i.e., the case of rational Weyl functions), and one can apply procedures, 
which are different from the procedures for the general-type case. 

Thus, we prove stability of solving inverse problems for matrix Dirac 
systems on semi-axis with both discrete and continuous components of spec¬ 
trum. Riccati equations play in essential role in the explicit solving of inverse 
problems, and so the classical works on Riccati equations by Leiba Rodman 
and coauthors are actively used in this paper. 

Explicit solutions of inverse spectral problems can be obtained either by 
applying procedures for general-type systems to some specihc spectral data 
(e.g., to rational Weyl or scattering functions) or by using several specihc for 
explicit solutions procedures. The hrst (general-type) approach was used, for 
instance, in [1,16,42] and [2, Sect. 6]. The second approach includes Crum- 
Krein method [8,29], commutation methods [9,17,18,27] and some versions 
of Backlund-Darboux transformation. Here we use our GBDT version of the 
Backlund-Darboux transformation (see [38,40,43] and references therein). 
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see also [4,5,33,48] and references therein on various versions of Backlund- 
Darboux transformations. 

In the next section. Preliminaries, we present some basic notions from 
system theory and formulate several results on Weyl functions. We also 
present GBDT procedure to solve inverse problem for systems (1.1) (more 
precisely, to recover self-adjoint Dirac systems from Weyl functions). Sec¬ 
tion 3 is dedicated to the proof of stability of this procedure. 

As usual, M stands for the real axis, C stands for the complex plane, C+ 
is the open upper half-plane { 2 : : A {z) > 0}, C_ is the open lower half-plane 
{z : ^(z) < 0 }, and the notation diagjdi,...} stands for the diagonal (or 
block diagonal) matrix with the entries di,... on the main diagonal. By ||A|| 
and by cr(A), we denote the /^-induced norm and the spectrum, respectively, 
of some matrix A. We say that the matrix X is positive (positive dehnite) 
and write X > 0 if W is Hermitian (i.e., X = X*) and all the eigenvalues of 
X are positive. 

2 Preliminaries 

2.1 Rational functions 

Recall that a rational matrix function is called strictly proper if it tends to 
zero at inhnity. It is well-known [25,30] that such an m 2 x mi matrix function 
ip can be represented in the form 

ip{z)=C{zh-A)-^B, (2.1) 

where .4. is a square matrix of some order n, and the matrices B and C are 
of sizes n X mi and m 2 x n, respectively. The representation (2.1) is called a 
realization of 99 , and the realization ( 2 . 1 ) is said to be minimal if n is minimal 
among all possible realizations of ip. This minimal n is called the McMillan 
degree of ip. The realization (2.1) of ip is minimal if and only if 

n—1 n—1 

span Im A’^B = C^, span Im {A*)’^C* = C”, n = ord(4,), ( 2 . 2 ) 

k=0 k=0 

where Im stands for image and ord(4.) stands for the order of A. If for a pair 
of matrices {4., B} the hrst equality in ( 2 . 2 ) holds, then the pair {A, B} is 
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called controllable. If the second equality in (2.2) is fulfilled, then the pair 
{C, A} is said to be observable. 

Now, let matrix functions (f be contractive, that is, let (p*(p < (or, 
equivalently, < 1 ^ 2 ) hold. From [30, Theorems 21.1.3, 21.2.1] (see also 
[13, p. 191]), the next proposition easily follows. 

Proposition 2.1 Assume that is a strictly proper rational matrix func¬ 
tion, which is contractive on M and has no poles in C+, and let realization 
(2.1) be its minimal realization. 

Then, there is a positive solution X > 0 of the Riccati equation 

XBB*X + i{A*X - XA) + C*C = 0. (2.3) 

Clearly, under conditions of Proposition 2.1, (p{z) is contractive on C+ U M. 

In the case of the skew-self-adjoint Dirac system, we obtain Ricatti equa¬ 
tion with minus before BB*: 

XC*CX + i{AX - XA*) - BB* = 0. (2.4) 

This case should be dealt with in a different way and we shall do it separately 
(in the next paper). 


2.2 System (1.1): Weyl function and inverse problem 

Recall that Y{x,z) is the normalized by Y{0,z) = Im fundamental solution 
of Dirac system (1.1), where j and V have the forms (1.2). 

Definition 2.2 An m 2 x mi matrix function p{z) {z G C+) such that 


[Im^ (p{z)*]Y{x,z)*Y{x,z) 


}p{z) 


dx < 00 


(2.5) 


is called a Weyl function of the Dirac system (1.1) on [0, cx)). 


Remark 2.3 Aeeording to [13, Sect. 2] and [43, Sect. 2.2], the Weyl fune- 
tion (p(z) of the Dime system (1.1) always exists and is unique. Moreover, 
ip{z) is holomorphic and contractive in C+. 
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If Lp is rational, it can be prolonged (from C+) on M and C_ in a natural way. 
Each potential v corresponding to a strictly proper rational Weyl function is 
generated by a fixed value n G N and by a quadruple of matrices, namely, 
by two n X n matrices a and 5*0 > 0 and by n x ruk matrices (/c = 1, 2) 
such that the matrix identity 

aSo - Soa* = i(didt - (2.6) 

holds. Such potentials v have the form 

v{x) = -2id*e^^“‘5(x)-^e^^“d2, (2.7) 

pX 

S{x) = So+ Ait)A{t)*dt, A{x) = e^^^ds] . (2.8) 

Jo 

Definition 2.4 [13, 19] The potentials v generated by the quadruples 

{q;, S'o, di, ^ 2 } {where Sq > 0 and (2.6) holds) via equalities (2.7) and (2.8), 
are called pseudo-exponential. 

Theorem 2.5 [13] Let Dirac system with a pseudo-exponential potential v 
be given on [0, cxd) and let v be generated by the quadruple {a, S'o, di, ^ 2 }. 
Then the Weyl function tp of this system has the form 

p,{z) = -\'dlS]]\zIn - 9 = a- idid^S'o"^ (2.9) 

The following theorem (i.e., [13, Theorem 3.4]) presents a procedure of ex¬ 
plicit solution of the inverse problem (see also [22, Theorem 5.2] for the 
mi = m 2 case), which is basic for this paper. 

Theorem 2.6 Let p{z) be a strictly proper rational matrix function, which 
is contractive on M and has no poles in C+. Assume that (2.1) is its minimal 
realization and that X > 0 is a solution of (2.3). 

Then p{z) is the Weyl function of the Dirac system (1.1), the potential 
V of which has the form (2.7), (2.8), where the quadruple {a, 5*0, di, ^ 2 } is 
given {in terms of A, B, C and X) by the relations 

a = A + \BB*X, So = X-\ =B, A 2 = -iX-^C*. (2.10) 
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In particular, the identity (2.6) easily follows from (2.3) and (2.10). The 
uniqueness of our explicit solution of the inverse problem is immediate from 
a much more general uniqueness result. 

Proposition 2.7 [41] The solution of the inverse problem to reeover system 
(1.1) from its Weyl function is unique in the class of Dirac systems with the 
locally square summable potentials. 

Remark 2.8 We note that there are many quadruples generating the same 
pseudo-exponential potential. The quadruples, which are recovered using (2.10), 
have an important additional property: controllability of the pair 
This property is immediate from the controllability of the pair {A, B}. 

Furthermore, the matrices A, B and C in the minimal realizations (2.1) of cp 
are unique up to basis (similarity) transformations: 

i:=r-Ur, c = cr, b = t~^b, (2.11) 

where T are invertible m x m matrices. Choosing the realization of p with 
A, B and C instead of A, B and C, and adding the sign in the notations 
of the corresponding matrices a, di and X, we derive 

a = T-^aT, di = T-^di (i = l,2), X = T*XT. (2.12) 

Setting T = X~^^‘^U*, where U is unitary, we have X = 1^. Hence, (2.6) 
takes the form a — a* = i('di'd^ — ' 112 ^ 2 ). Moreover, for the case mi = m 2 = p, 
it was shown in [21] that U may be chosen in such a way that we have the 
block representations: 

/3 :=a- idi(di + ^ 2 )* = 

where 

C = C* = diag{fi4,,t242,---,44j, a0)eC_, (2.14) 

and u is some n x p matrix, n := 


13 0 

0 C 


A = 


A 

u 


^2 = 


^2 

—U 


(2.13) 


6 



Now, introduce Dirac operator H associated with the differential expres¬ 
sion 

x) = /(^)’ ( 2 - 15 ) 

the domain of which consists of all absolutely continuous functions / from 
L^(0, cxd), such that Tif G -h^(0, cxd) and the initial condition 

[/, -4] /(O) = 0 

holds. Using (2.13), it is shown in [21] (see also [22, Sect. 2]) that the real 
eigenvalues of H are concentrated at the points 4 and have multiplicities Uk, 
whereas the continuous spectrum of H is described by /3, and ^ 2 . Namely, 
the spectral density g of H has the form 

g{t) = g{t)*g{t), g{t) := Ip - i(di ATitIn-n - (2.16) 

In view of the mentioned above connections between the quadruple 
{a, 5*0, di,'^ 2 } and the corresponding Weyl and spectral functions, we can 
consider this quadruple as some spectral data. 

3 Stability of explicit solutions 

1. The following lemma is a stronger statement than Proposition 2.1. (The¬ 
orem 7.4.2 from [30] is used for its proof in addition to the Theorems 21.1.3 
and 21.2.1 from [30], which yield Proposition 2.1.) 

Lemma 3.1 [15] Assume that a strictly proper rational m 2 x mi matrix 

function (f{z) is contractive on M, and that (2.1) is its minimal realization. 

Then there is a unique Hermitian solution X of the Riccati equation (2.3) 
such that the relation 



a{A + iBB*X) (3.1) 

holds. This solution X is always invertible. It is also positive if and only if 
ip{z) is contractive in C+. 
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Further, in our procedure to recover the potential v from the Weyl function 
99 , we shall look for this particular solution X of (2.3). More precisely, we 
start with the strictly proper rational m 2 x mi matrix function ip{z), which 
is contractive on R and has no poles in C+. Hence, ip{z) is contractive in C+, 
and so, according to Lemma 3.1, we have X > 0. By Qn we denote the class of 
triples {A, B,C} which determine minimal realizations (p{z) = C{zln — A)~^B 
of m 2 X mi matrix functions <f{z) contractive on R U C_|_. First, we consider 
the stability of recovery of X from {A, B,C} G Qn- 

Definition 3.2 The recovery of X >0, satisfying (3.1), from the minimal 
realization (2.1) {where {A,B,C} G Qn) of i^{z) is called stable if for any 
e > 0 there is S > 0 such that for each {A, B,C} satisfying conditions 

{A,B,C}egn, \\A-A\\ + \\B-B\\ + \\C-C\\ <S (3.2) 

there is a solution X = X* of the equation XBB*X + \{A*X — XA) +C*C = 0 
in the neighborhood ||X — X|| < e of X. 

The stability of solutions X of an important class of Riccati equations was 
shown in [36, Theorem 4.4] for a somewhat wider class of perturbations than 
described in our dehnition and we shall use this theorem in order to prove 
our hrst stability statement. 

Theorem 3.3 The recovery of X >0, satisfying (3.1), from the minimal 
realization (2.1) of i^{z) {with {A,B,C} G Qn) is stable. 

Proof. Assuming that a minimal realization (2.1) of ‘f{z) is given (that is, 
matrices A, B and C are given), we consider equation (2.3). Putting 

Ho = -i{A + cin) (c G R), Co = -C*C, Do = BB*, (3.3) 

we see that the equation (2.3) coincides with the Riccati equation XDqX + 
XAo + AqX — Co = 0 considered in [36, Subsection 4.2]. 

Now, we deal with the conditions (i)-(iv) (on the coefficients Ao, Co, Do) 
from [36, Subsection 4.2]. (Only perturbations satisfying these conditions 
are allowed in [36, Subsection 4.2] and we will show that the conditions 
(i)-(iv) are satisfied in the case {A,B,C} G Qn-) Equalities (3.3) and the 



fact that the pair {A,B} is controllable imply that conditions (i) and (ii) 
in [36, Subsection 4.2] are fulhlled. In a similar way we derive that conditions 
(i) and (ii) are fulhlled for the Ricatti equations XDqX+XAq+AIX — Cq = 0 
corresponding to all the triples {A,B,C} G Qn- For sufficiently large values 
of |c|, the requirement (iii) that the matrix 


-Co RS 

^0 Dq 


(3.4) 


satishes the condition det H ^ 0 and that signature of H equals zero is also 
fulhlled. Clearly, c may be chosen so that (iii) is valid if we substitute the 
triple {A,B,C} with any triple {A,B,C} G Qn in some small neighborhood 
of {A,B,C}. Finally, according to Lemma 3.1, there are hermitian solutions 
of equations XDqX + XAq + A^X — Cq = 0, that is, condition (iv) holds. 
Since conditions (i)-(iv) from [36, Subsection 4.2] are fulhlled, the stability 
with respect to perturbations in the class Qn will follow from the stability in 
the sense of [36, Theorem 4.4]. 

Using again Lemma 3.1, we choose the solution X > 0 of (2.3) satisfying 
(3.1). It is immediate that one of the equivalent statements from [36, Theo¬ 
rem 4.4] is valid for our X. That is, according to (3.1) and (3.3), the equality 
Q'(A) = 0 holds for each A from the set 


a(i(^ + DoX)) n a( - i(A* + XD^)) 

= a{A + iBB*X + cln) n a{{A+iBB*X + cln)*), (3.5) 


and so the statement (d) from [36, Theorem 4.4] holds. Therefore, by virtue 
of [36, Theorem 4.4], X is a stable and isolated solution of (2.3). ■ 

2 . Now, we will show that small perturbations of the quadruple {a, So, '^ 1 ,^ 2 } 
result in small perturbations of the corresponding potential v. We note that 
we consider only perturbations which do not change mi, m 2 and n . 


Definition 3.4 The quadruple {a, S'o, di, ^ 2 } is ealled admissible i/S'o > 0 
and (2.6) holds, and it is called spectral if it is admissible, the pair {Q:,'di} 
is controllable and 


a{a) C (M U C_). 


(3.6) 
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Remark 3.5 Theorem 2.5 and Remark 2.3 show that the Weyl funetion cor¬ 
responding to any pseudo-exponential potential is rational and contractive. 
Then Theorem 2.6, Proposition 2 .7 and Lemma 3.1 imply that this potential 
{uniquely recovered from the Weyl function) is generated, in particular, by 
a spectral quadruple. In other words, each pseudo-exponential potential is 
generated by some spectral quadruple. 

Theorem 3.6 Let a spectral quadruple {a, S'o, 't? 2 } be given. Then, for 
any £ > 0 there is 5 > such that each pseudo-exponential potential v 
generated by an admissible quadruple {5, S'o, 'f? 2 } satisfying condition 

||q! - 511 + \\So - S'oll + ll^?! - + ||^?2 - 'f? 2 || < d 

belongs to the e-neighborhood of v generated by {a, S'o, 't? 2 }, that is, 

sup ||u(a;) — u(a;)|| < £. (3.7) 

x€[ 0 ,oo) 

In order to prove the theorem above we generalize (for the case when nii does 
not necessarily equal m 2 and 5*0 does not necessarily equal In) some results 
from [ 22 ] on asymptotics of 

pX 

Q{x):=So + 2 (3.8) 

Jo 

Lemma 3.7 The following relations are valid for a spectral quadruple 

{a,S'o,?9i,^?2} : 


lim Q{x)-^ = 0, lim ||Q(x)-= 0. (3.9) 

X—^OO X—^OO 

Proof. The proof uses some steps from the proof of [22, Proposition 3.3]. 
Since Q{x) is increasing and is bounded below by 5*0 > 0, there is a limit 
kq := lima;_j.oo Q(x)“^. Next, we prove that kq = 0. From the dehnition 
(3.8) and identity (2.6) we derive 

aQ{x) — Q{x)a* = aSo — Soa* — i (e^^'^“d 2 'd 26 ~^^'^“* — ^? 2 'd 2 ) 

= idWl - ie2"“^?2^?2e“""“‘. (3.10) 
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Multiplying (from both sides) the left-hand side and right-hand side of (3.10) 
by we obtain 

Q{x)~^a — a*Q{x)~^ — iQ{x)~^'&i'&lQ{x)~^ 

= -iQ(x)-^e2^^"d2^9^e-2^^“*g(x)-h (3.11) 

Passing in (3.11) to the limit, we see that 

lim g(x)“^e^'^"?92^?2e"^‘'^“*g(a^)"^ = Ki^qa - a*KQ - ikq'&i'&Ikq). (3.12) 

X^OO 

On the other hand, formula (3.8) yields 

^Q(x)-^ = -2Q(x)-^e2'^“d2d;e-2““*g(x)-\ 
ax 

and so we have 

g(f)"^e^'‘“d2^?2e“^'‘“*g(f)"^df = - kq) < oo. (3.13) 

Taking into account (3.13) and the fact that there exists a limit of the expres¬ 
sion integrated in (3.13) (see (3.12)), we derive that this limit equals zero. 
That is, we rewrite (3.12) in the form 



Kqa - a*KQ-\HQ'di'&\KQ = (3.14) 

Moreover, since the left-hand side in (3.12) tends to zero, the second equality 
in (3.9) is already proved. 

Recall that the first equality in (3.9) is equivalent to kq = 0. Now, we 
prove kq = 0 by negation. For this, we rewrite (3.14) in the form a*KQ = 
— which implies that the range of kq is an invariant subspace 

of a*. Thus, assuming kq ^ 0, we obtain that there is an eigenvector kqq of 

a*\ a*KQg = CKqg, Kqg ^ 0, g E C”. 

Finally, consider the expression ig*{KQa — a*KQ)g. In view of (3.6), for 
the eigenvalue c of a* we have ^^(c) > 0, and so 

ig*{KQa - a*KQ)g = i(c - c)g*HQg > 0. 
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On the other hand, we have '&*iKQg ^ 0 because the pair {ajdi} is control¬ 
lable. Hence, the inequality \g*{KQa — a*KQ)g = —g*KQ'di'dlKQg < 0 follows 
from (3.14). We arrive at a contradiction, that is, kq = 0. ■ 

In the case of admissible quadruples {a, 5*0, di, ^ 2 }, the matrix identity 

aS{x) — S{x)a* = iA(x)jA(a;)* (3.15) 

(see [13, formula (3.6)]) coincides with (2.6) at a; = 0 and easily follows from 
(2.6) and (2.8) for a; > 0. In other words, a, S{x) and A(x) form an 5'-node 
(and, moreover, the so called Darboux matrix function corresponding to v{x) 
coincides with the transfer matrix function [43-45] in Lev Sakhnovich sense). 
Using (2.8), (3.8) and (3.15), we derive 

0'(i) = (e‘“S(i)e-‘“-)', Q(0)=S(0) ([q':= , 

and so the following equality is valid: 

Q(x) =e^^“5(x)e-'^“‘. (3.16) 

Proof of Theorem 3.6. Now, we consider a pseudo-exponential potential v 
generated by the spectral quadruple {a, S'o, di,'d 2 } and pseudo-exponential 
potentials v generated by admissible quadruples {5, /Sq, "di,-(92} belonging to 
a neighborhood of {a, S'o, di,'d 2 }. The matrix function Q corresponding to 
{5, ^o, di, d 2 } is denoted by Q. In view of (2.7) and (3.16), we have: 

v{x) = -2id*Q(a;)-^e2'""d2, T{x) = -2id*Q(a;)-^e2“"d2. (3.17) 

It is immediate from the proof of Lemma 3.7 that (3.11) holds for admissible 
quadruples as well. That is, we may rewrite (3.11) for {5, 5*0, di, d 2 }: 

Q(x)~^a — a*Q{x)~^ — iQ{x)~^'6i'&lQ{x)~^ 

= -ig(a;)-^e2^^"d2d;e-2^^"*g(x)-h (3.18) 

Since Q and Q are monotonic and the hrst relation in (3.9) is valid, we may 

choose xo > 0 and some neighborhood of {a, 5*0, di, d 2 } so that Q{x) and 
Q{x) are large enough for x > Xq, and so the left-hand sides of (3.11) and 
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(3.18) are small enough. Hence, the right-hand sides of (3.11) and (3.18) are 
also small enough. Therefore, taking into account (3.17), we see that for any 
e > 0 there are Xq > 0 and > 0 such that the next inequality holds in 
the 5i-neighborhood (||q; — 5|| -|- Us'd — S'o|| -|- ||di — di|| -|- ||'d 2 — d 2 || < (^i) of 
{a, S'o,di,d 2 }: 

sup ||x(x) — h(x)|| < e. (3.19) 

a:G[xo,oo) 

It easily follows from the dehnition of Q and Q and from (3.17) that there is 
some 52-neighborhood of {a, S'o, di, ^ 2 }, where we have 

sup \\v{x)— v{x)\\ < e. (3.20) 

a;€[0,xo) 

Clearly, inequalities (3.19) and (3.20) yield (3.7) (for 6 = min(5i,52)). ■ 

Remark 3.8 It is immediate from the second relation in (3.9), the first re¬ 
lation in (3.17) and Remark 3.5 that all pseudo-exponential potentials tend 
to zero at infinity. 

3. Theorems 2.6, 3.3 and 3.6 as well as Lemma 3.1 yield the result below on 
the stability of the procedure of solving inverse problem. 

Theorem 3.9 The procedure {given in Theorem 2.6) to uniquely recover the 
pseudo-exponential potential v of Dirac system (1.1) from a minimal real¬ 
ization of the Weyl function (he., of some strictly proper rational m 2 x mi 
matrix function, which is contractive in C+) is stable once we agree to choose 
such a positive solution X of the Riccati equation (2.3) that (3.1) holds. 
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